Rapid and robust spin state amplification 
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Electron and nuclear spins have been employed in many of the early demonstrations of quantum technology 
(QT). However applications in real world QT are limited by the difficulty of measuring single spins. Here we 
show that it is possible to rapidly and robustly amplify a spin state using a lattice of ancillary spins. The model 
we employ corresponds to an extremely simple experimental system: a homogenous Ising-coupled spin lattice in 
one, two or three dimensions, driven by a continuous microwave field. We establish that the process can operate 
at finite temperature (imperfect initial polarisation) and under the effects of various forms of decoherence. 

PACS numbers: 



The standard approach to implementing a quantum technol- 
ogy is to identify a physical system that can represent a qubit: 
it must exhibit two (or more) stable states, it should be ma- 
nipulable through external fields and possess a long decoher- 
ence time. Provided that the system can controllably interact 
with other such systems, then it may be a strong candidate. 
Electron and nuclear spins, within suitable molecules or solid 
state structures, can meet these requirements. However the 
drawback with spin qubits is that they have not been directly 
measured through a detection of the magnetic field they pro- 
duce. The magnetic moment of a single electron spin is orders 
of magnitude too weak to be detected by standard ESR tech- 
niques and even the most sensitive magnetometers still fall 
short of single spin measurement [1] - meanwhile the situation 
with nuclear spins is worse still. In a few special systems it 
is possible to convert the spin information into another degree 
of freedom. For example, a spin-dependent optical transition 
allow spin to photon conversion in some crystal defects 12- 
4 1, self-assembled semiconductor quantum dots J5] |6), and 
trapped atoms held in a vacuum Q. Alternatively, spin to 
charge conversion is an established technology in lithographic 
quantum dots |8 1. However, the majority of otherwise promis- 
ing spin systems do not have such a convenient property [9| 
and therefore cannot be measured directly. 

One suggested solution is to 'amplify' a single spin, by us- 
ing a set of ancillary spins that are (ideally) initialised to |0). 
We would look for a transformation of the form 

|0) \of n -> |0) \of n |1) \of n ->■ |1) \lf n , (1) 

the idea being that the n ancillary spins constitute a large 
enough set that state of the art magnetic field sensing tech- 
nologies can detect them. Note that the transformation need 
not be unitary or indeed even coherent: since the intention is 
to make a measurement of the primary spin, it is not necessary 
to preserve any superposition (that is, we need not limit our- 
selves to transformations that take a |0) \$f n + ft |1) \0)® n 
to a cat state like a |0)® n+1 + /3 |l)® n+1 ). 
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This is a rather broadly defined transformation and there 
are a number of ways that one might perform it. Clearly 
one would like to find the method that is the least demanding 
experimentally. Previous authors have proposed schemes us- 
ing a strictly one-dimensional (ID) homogeneous lattice with 
continuous global driving |10|, and an inhomogeneous three- 
dimensional (3D) lattice with alternating timed EM pulses 
ITTl . The former result has the advantage of simplicity but the 
rate at which amplification occurs will inevitably be limited 
by the single dimension of the array; moreover such a sys- 
tem must be highly vulnerable to imperfect initialisation (i.e. 
finite temperature). Here we generalise to a homogeneous 
two-dimensional (2D) square lattice, showing that a contin- 
uous global EM field can drive an amplification process that 
succeeds at finite temperatures (imperfect initialisation of the 
ancilla spins) and in the presence of decoherence. By bringing 
the global EM field onto resonance with certain transitions, we 
are able to create a set of rules that govern locally how spins 
propagate over the lattice. We then look at the rate of increase 
in the total number of flipped spins as a measure of quality of 
the scheme. While our focus is on the 2D case, we are also 
able to predict the performance of the amplification protocol 
for a homogeneous 3D lattice with continuous driving. 

The case of a ID lattice has been studied in detail by Lee 
and Khitrin 1 10 1. Before moving to the 2D spin lattice that will 
form the core of the paper, we first recall how to simplify the 
description of this (semi-infinite) ID spin chain, with nearest 
neighbour Ising (ZZ) interactions. Under a microwave driving 
field of frequency to, the Hamiltonian is given by 

oo 

V. = e t a l z + .halaf 1 + 2a t a l x cos(wt) (2) 

i=l 

€i is the on-site Zeeman energy of spin i, and Ji is the magni- 
tude of the coupling between spins i and i + 1. £1 describes the 
coupling of spin i and the microwave field. In this case, spin 
i = 1 is the one whose state is supposed to be amplified. If 
we assume that the chain is uniform, such that fij = f2, e; = e 
and Ji = J, then moving into a frame rotating at frequency 
u), making a rotating wave approximation and setting lu = e 
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leads to 

oo 

H = J*>i +1 + *K- (3) 

i=l 

In order to understand the dynamics of the system, is it in- 
structive to explicitly separate all terms that involve a particu- 
lar spin k: 

H = J{cj k z - 1 +a k + 1 )a k z +na k x + ^ tta x +Ja l z cj l z +1 +tt(T k 

i^{fc,fc-i} 

(4) 

Choosing a driving field such that £1 <C J means that spin k 
will only undergo resonant oscillations when the first term in 
Eq.|4]goes to zero - i.e. when the two spins neighbouring spin 
k are oriented in opposite directions. In any other configura- 
tion the Ising coupling takes the spin k off resonance with the 
microwave and no appreciable dynamics are expected. 

Let us now define a subset of states S that exist in the spin 
chain Hilbert space, \n), which have the first n spins of the 
chain in state \\) with the rest \\). If the rule we just derived 
holds exactly these states define a closed subspace. We may 
then write a very simple isolated Hamiltonian for this sub- 
space: 

oo 

Hs = nJ2\n)(n + l\ + \n+l)(n\. (5) 

71=1 

With this simplification of the ID Hamiltonian in mind, 
we progress now to a semi-infinite square spin lattice with 
nearest-neighbour ZZ interactions. For this case we have 

oo oo 

U = Y^£u e(ri' j +Jai' j ai +1 ' j +Jai' j ai' j+1 +2na^ cos{u>t). 

i=l 1 = 1 

(6) 

By again considering the terms affecting a particular spin in 
the main body of the lattice (fc(> 1), l(> 1) say) we find for 
lu = e and after moving to a rotating frame and making the 
rotating wave approximation: 

n = Ja k ' l (a k+1 ' 1 + a k z l+1 + a k ~^ 1 + a k z ' l ~ l ) + ... (7) 

where we do not explicitly write out terms not involving spin 
(k, I). The microwave is now only resonant for spin (k, I) if it 
has two neighbour spins in each orientation. For a spin on the 
edge of the lattice there are an odd number of neighbours so 
resonance cannot be achieved in this case. However, applying 
a second microwave with oj = e — J allows resonant flips on 
the edge if two neighbours are down and one up - and this 
second field has no effect on the bulk spins. 

The spin to be measured is the corner spin (i = j = 1) 
and so would form part of a wider computational apparatus. 
We may therefore assume that it is a different species with a 
unique resonant frequency. The dynamics of the whole lattice 
may then be summarised by three rules (in order of prece- 
dence): 

1. The corner (test) spin is fixed. 




FIG. 1: Partition states arranged into a lattice. Edges represent a 
coupling through the Hamiltion of strength fi. Weights represent the 
number of different paths through the lattice to a given state. 

2. An edge spin can flip if it has one of its neighbours up 
and two down. 

3. A body spin can flip if it has two of its neighbours up 
and two down. 

We begin by supposing all spins are initialised in the 'down' 
state apart from the test spin, which is located in the upper left 
hand corner of our lattice. We can describe this initial state by 
choosing two basis elements: |0) when the test spin is down, 
and |1) when the test spin is up. Using our heuristic rules we 
can see that these two states do not couple to each other - that 
(0| H |1) = 0. In fact |0) does not couple to any other state, so 
if we start in the |0) state no amplification occurs, as desired. 

We will now seek to construct a basis for the subspace con- 
taining our system evolution, by looking at states connected 
by our Hamiltonian. It will be convenient to represent these 
states on the nodes of a graph, using the edges to represent 
non-zero elements of the Hamiltonian. 

Our starting point is the state |1), with just the comer spin 
'up'. From this position our rules allow two possibilities: ei- 
ther the spin to the right of the corner flips, or the spin below 
it flips (see Fig. [T}. In each case the magnitude of the transi- 
tion matrix element is O. As we continue this procedure, we 
notice that the states that arise for each excitation number can 
be characterised by a non-increasing sequence of integers that 
represent the number of 'up' -spins in each column of the lat- 
tice (see Fig. [T]). Such sequences can also be used to define 
partitions of at integer: ways of splitting an integer up into a 
sum of other integers, e.g. 3 = 3 = 2 + 1 = 1 + 1 + 1. 
In fact, the states that arise are in 1-to-l correspondence with 
such partitions; we call these states 'partition states' and de- 
note them with standard partition notation (see Fig. [TJ. The 
graph we have just described is depicted in Fig. [T] is known 
as 'Young's lattice' and arises in areas of pure mathematics, 
such as the representation theory of the symmetric group, and 
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the theory of differential posets. We have drawn weights be- 
neath each state, recording the number of ways the state can 
be constructed. We will now further reduce the dimension of 
this basis by eliminating combinations of states which are in- 
accessible. 

Starting with |1) we see that 

(l|fl"(ai,i|VM)+«al^a>) = tt(ai,i+a 2 ) so |1) 
does not couple to the two-excitation state — IV^)- We 

can eliminate this, leaving a single orthogonal, coupled state 
with two excitations: \2) :— dV'i.i) + IV^))- 

We may continue to build up coupled states with larger ex- 
citation numbers, and in fact we find that there is only a sin- 
gle coupled state in each case (i.e. we can always eliminate 
k — 1 combinations of partition states with k excitations). To 
see this, first suppose we have the coupled state with k exci- 
tations, which by analogy with the ID case we write as |fc). 
We can write \k) = ^ Sieprfc) c « IV , *)> where P(k) is the 
set of partitions of the integer k and Nj. a normalisation fac- 
tor. We want to construct the state \k + 1) by eliminating the 
fc-dimensional subspace with k + 1 excitations, to which 
does not couple. 

Let = J2jeP{k+i) a j IV'j) anc l consider the states \tp) 
such that 

= {k\H\4>)= ]T E 4 a AHHWj) 
ieP(k) jeP(k+i) 

but (ipi \H\ipj) = SI if z is a parent of j (a state connect to j, 
in the lattice row above it), and otherwise, so 

= {k\H\t()) = a i 
jeP(k+i) 

This is the equation of a hyperplane in \P(k + 1)| dimen- 
sions, defining the states that are not coupled to \k) through 
the Hamiltonian. There is a unique single state orthogonal to 
this hyperplane, fa = £ iepore „ ts(i) Cj, to which \k) couples. 
So the only state with k + 1 'up'-spins that couples has 
coefficients proportional to (3j. After normalisation, we call 
this state \k + 1). 

Unfortunately there is no easy way to write down the par- 
tition states and weights for the nth row of the lattice. Fortu- 
nately, for our purposes, we only need to know that the states 
|fc) exist and what the coupling between them is. To find this 
coupling, consider 
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To find the N n we need the sum of the squares of the weights 
of partitions in a given row. A standard result about Young's 
lattice immediately gives us this sum: n\ lfl2l . 
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FIG. 2: Expected total polarisation against time. Time in units of 
dephasing rate r = 1. The gradient of the 'one dimension with 
decoherence' line tends to | asymptotically. 



Referring back to Eq. (|8j), and using Ni 

H = Vn(\n- 1) (n| + \n) 



we see that 



(n-l| 



(9) 



In essence we have established a linear sequence of states, 
each coupled to the the next analogously to the states on a ID 
chain [5] However, each of our states is in fact a superposi- 
tion of many configurations of the 2D array, and crucially the 
effective coupling from each state to the next increases along 
the sequence. 

It has been shown (e.g. [1 1 3 '] ) that a quantum state released 
at the end of a semi-inifinite chain of states, with constant 
couplings, will travel ballistically: the average position of the 
state along the chain is proportional to the time passed, and 
inversely proportional to the coupling strength. Since, in the 
one-dimensional case, the position is proportional to the num- 
ber of spins that have flipped, we have that the total polarisa- 
tion will increase linearly with time. 

We can establish the rate of propagation in the 2D case 
using the ansatz that the time taken to travel between two 
neighbouring nodes is inversely proportional to the strength 
of the coupling between them. The total time is then t^D oc 
SILi ~^f\ — « 3 ■ As in the one-dimensional case, the position 
along the chain corresponds to the the number of spins that 
have flipped, and so we would expect the total polarisation to 
be proportional to t 2 . This prediction of a quadratic speed- 
up of signal going from ID to 2D is the central result of our 
paper, and was confirmed by simple numerical simulations of 

Unfortunately the mapping from 2D to ID is not readily 
extendible to 3D. However, our results so far could have been 
anticipated using simple dimensional arguments; if one pos- 
tulates that the rate of spin propagation is proportional to the 
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boundary of the region, one can predict the correct scaling be- 
haviour. In ID the boundary size is independent of the region 
size; no matter how many spins have flipped, it still has size 
one. The coupling strength between states \n) is constant. In 
the 2D case, the boundary size scales with the square root of 
the area, and the coupling goes with *fn. In 3D, the bound- 
ary scales like the cube root of the volume squared, and so we 
expect the coupling to scale as ni . Following similar logic to 
that used in 2D case: t^n oc y^™_i \ — , and so n ~ i 3 . 

We now consider the effect of decoherence. Much of 
the early work on continuous time quantum random walks 
looked at the speedup they afforded over their classical coun- 
terparts IT41 . but didn't make any statement about the con- 
ditions under which we would expect the quantum walk to 
exhibit classical behaviour, as we might expect in a regime of 
suitably heavy dephasing, say. 

We begin by considering a collective noise operator: L = 
T n n \ n ) ( n \- This represents noise that applies uniformly to 
the whole lattice: global fluctuations in the magnetic field, for 
example. As the effect of this type of noise depends only on 
the number of 'up' spins, the system remains in the reduced 
basis of number states calculated earlier, with only the coher- 
ences between these states affected. 

Our starting point is the Lindblad master equation 

p = i[p,H} + ~r (2LplS - iSLp - pL^L) . (10) 

We proceed by splitting up the equation into diagonal and off- 
diagonal terms: 
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(12) 



where g^ is the coupling between states i and j. In the limit 
of heavy dephasing (r ^> g), we have a process similar to 
adiabatic following, and we can make the approximation 



^ Pikdk] - ^ Pkj9ik 



V *=±J 



fe=±l 



n(n— 1) 



variables and solve for 



We consider the p^j as a set of , 
them in terms of the pa. Neglecting terms that are second 
order in and substituting back into Eq. (Hi gives 



Pa = 



E 



2|flijl 



{Pii ~ Pjj) 



Our quantum chain formally reduces to a classical Markov 
chain on the same statespace, with transition rates propor- 
tional to the coupling squared. 

Although states with more 'up' spins decohere more 
quickly, the decoherence rate T is not multiplied for higher 



states, as it is the relative decoherence rate between neigh- 
bouring states, which is of importance. 

In one-dimension gry = 1 and we are reduced to a simple 
random walk on a semi-infinite line. By analogy with simple 
diffusion we expect that the resulting distribution is roughly 
Gaussian, with the expected number of flipped spins going 
with \fi: the rate of spin propagation drops from t to \fi. This 
result was confirmed numerically (Fig. |2j. 

In the two-dimensional case g^ = = i + 1: We get 

a random walk with increasing transition rates. Numerically 
(Fig. [2}, we find that the rate of spin propagation drops from 
t 2 to t - still an encouraging scaling. 

We also investigated the 'individual noise' case, where the 
dephasing occurs independently on each site (see supplemen- 
tary material). Although the calculations differ, we see the 
same rate behaviour as for the collective noise case. 

Finally we consider imperfect inital polarisation (i.e. finite 
temperature) - a property exhibited by any real experimen- 
tal system. As discussed in the supplementary material IPT31 . 
a fortuitous consequence of the propagation rules is that our 
system is particularly robust against this source of error; be- 
low an initialization threshold of approximately 4% it is ex- 
tremely unlikely that a false positive will occur. This places 
our protocol well within experimental capabilities; for exam- 
ple for an array placed in a standard W-band electron spin 
resonance system (100 GHz) and cooled using liquid 4He to 
1.4 degrees Kelvin, only 3.1% of electron spins will be in the 
'up' state. 
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